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Abstract 

The spread between two lines in rational trigonometry replaces the concept of angle, allowing the complete 
specification of many geometrical and dynamical situations which have traditionally been viewed approximately. 
This paper investigates the case of powers of a rational spread rotation, and in particular, a curious periodicity 
in the prime power decomposition of the associated values of the spread polynomials, which are the analogs in 
rational trigonometry of the Chebyshev polynomials of the first kind. Rational trigonometry over finite fields 
plays a role, together with non-Euclidean geometries. 

1 Introduction 

This paper investigates the role of the spread polynomials S n (s) of rational trigonometry in understanding a simple 
dynamical system arising in elementary geometry. Spread polynomials were introduced in [5J, and arise naturally 
when Euclidean geometry is studied algebraically, with angles replaced by spreads, and they play also a role in 
geometry over finite fields, and in fact over general fields not of characteristic two. Furthermore their importance 
extends to elliptic and hyperbolic geometries, see [5]. Spread polynomials also have interesting number theoretical 
properties that set them apart from Chebyshev and other orthogonal polynomials, see [6J and [2]. 

In this paper we study a phenomenon which occurs when we analyze iterates, or powers, of a particular rotation 
p in the Euclidean plane. If p is the rotation by an angle 9 which is a rational multiple of 7r, then the subsequent 
powers p n are easily understood. However when 9 is an irrational multiple of 7r, the situation exhibits chaotic 
features, and illustrates the 'butterfly effect': the future state of the system is determined by the accuracy by which 
we know the real numbers which specify initial conditions, and even a small error will eventually result in complete 
uncertainty — this is the standard view. 

With rational trigonometry, increasing uncertainty is replaced with increasing complexity. By measuring ro- 
tations p with a spread s, not an angle 9, the spread polynomials arise naturally from iterates, and the chaos 
disappears, replaced instead by an increasing escalation in the size and complexity of the numbers S n (s) that 
describe the evolution of the system. These numbers are far from random, and they tend to incorporate interesting, 
and indeed mysterious, number theoretic information. 

For a rational spread s, we will prove a theorem about the spreads S n (s) associated to the powers p n of the 
rotation p with initial spread s. Such rational spreads are common in geometry, for example any configuration 
of lines with rational equations has only rational spreads. Of particular interest is that we will need to examine 
finite geometries in order to understand the situation over the rational numbers, and that forms of non-Euclidean 
geometry also arise. Our work can also be viewed in the context of investigations of Chebyshev polynomials over a 
finite field, initiated already in [5]; see also [1] and [3]. 



1.1 Powers of a rotation 

To provide some motivation for the spread polynomials, recall that with the ISO size standards, an AA sheet of paper 
has proportions in the ratio \/2 : 1. So if you fold a piece in two lengthwise, the result has the same proportions 
but only one half of the area, and is called A5 size. The angle in degrees between the long side Iq and the diagonal 
h of an A4 sheet may be described by the formula 

9 ISO = f 180 arcsin vT/s) /tt. 
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This has a numerical value in degrees of approximately Oiso ~ 35.26. 

Suppose we rotate l\ by Oiso to get I2, then rotate I2 by Oiso to get I3, and so on. We could also say that l 2 is 
the reflection of Iq in l-y, that I3 is the reflection of l\ in I2, and so on. This gives a sequence li,h,- ■ • of concurrent 
lines whose respective angles with l form the sequence Oiso, 29iSOi ^iso-, ' ' ' • Since #/so is n °t a rational multiple 
of 7r, these sequences of lines and angles will never repeat, and so the lines l n are all distinct, as are the angles 
n0 iso . 




Figure 1: Multiples of Oiso, or equivalently s = 1/3 

Figure [T] shows the first sixteen such lines, exhibiting the famous three gaps phenomenon, which was originally 
a conjecture of Steinhaus, (see for example [3]): for any n, there are at most three different angles found between 
any two adjacent lines in the set {Iq, h, h, • • • In}- 

The angle between Iq and Iiq is approximately 

m ISO « 16 X 35.26 mod 360 = 564. 16 - 360 = 204. 16. 
We have lost information, however, and with the better initial approximation Oiso ~ 35. 264 389 682 we get instead 

me ISO « 204.230 235. 

But no matter what accuracy we have for Oiso, computing larger and larger multiples of it, mod 360, inevitably 
results in increasing error — until after a finite number of multiples all knowledge of the position is lost. This is a 
simple example of the well-known butterfly effect, a feature of a wide range of dynamical systems, and an inevitable 
consequence of working in the framework of classical trigonometry, which generally deals only with approximations 
to real values. The usual idea is that our understanding of the future evolution of many completely specified systems 
is limited by the precision with which we know the real numbers that appear in the initial conditions. 

With rational trigonometry, introduced in [BJ, we set our sights higher — we aim to describe such a system 
completely precisely, until we run out of computing power, memory space, or patience. The price we pay for 
more accuracy is — more complexity. For dynamical systems which can be expressed using rational numbers and 
polynomial transformations, the future evolution is completely known, but will generally be increasingly difficult and 
expensive to write down as time goes on. We will show that rotation by the angle Oiso exhibits this phenomenon, 
and brings out interesting number theoretical questions. 

It will be important for us to consider not just rotations in the Euclidean plane, but also over finite prime fields. 
For this it is often more convenient to work not with a unit circle, but rather with the associated projective line. 
The metrical geometry of the projective line is described by a projective version of rational trigonometry, which 
also extends to higher projective spaces, and it is in terms of this that we formulate rotations and reflections. 
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1.2 Basic rational trigonometry 



Here is a quick review of some main ideas from affine rational trigonometry in the plane; the main reference is [6], 
see also |10j . We work over a general field, not of characteristic two. The primary measurement is the quadrance 
between two points A\ = [xi,yi] and A 2 = [x 2 ,y 2 ], defined by 

Q(A U A 2 ) = (x a - xxf + (y 2 - yif . 

Over the real numbers, quadrance is the square of distance — or better yet distance is the square root of quadrance. 
However in rational trigonometry we wish to avoid square roots, and so distance plays no role. This allows the 
theory to extend to more general fields. 

The short side, long side and diagonal of an A4 sheet of paper form a right triangle with quadrances Q%, Q 2 and 
Q 3 in the ratio 1:2:3, and Pythagoras' theorem takes the form 

Qi + Q 2 = Q-3- 

In rational trigonometry the separation of two lines is measured by a spread, not an angle. A line I with equation 
ax + by + c = is a null line precisely when a 2 + b 2 = 0. The spread between two non-null lines with equations 

a%x + biy + c\ = and a 2 x + b 2 y + c 2 = 

is defined to be the number 

= («i&2 - ggji) 
S " (a? + 6?) («! + &!)• 

This formula also gives the spread between lines with direction vectors (ai, &i) and (02, b 2 ). Over the real numbers 
the spread s between two lines lies between and 1, being when lines arc parallel, 1 when lines arc perpendicular, 
and in general the square of the sine of an angle 9 between them. There are many such possible angles (e.g. 6, 
7r — 6, 2ir + 6, 3tt — 8, ■ ■ ■ ), but the square of the sine — the spread — is the same for all. 

The spread between two non-null lines can be expressed as a ratio of two quadrances — an opposite quadrance 
to a hypotenuse quadrance in a right triangle formed from those lines, so for example the spread between the long 
side of a piece of AA paper and its diagonal l\ is s(Zq,Zi) = s = 1/3. Unlike the notion of angle, spread really 
is defined between lines, not rays. Using spreads instead of angles makes much of the study of triangles — that is, 
trigonometry dramatically simpler. This is explained at length in [Sj. 

Define a number s to be a spread number precisely when s (1 — s) is a square in the field. The Spread number 
theorem ([5J Chapter 6]) asserts that a number s is the spread between two lines precisely when s is a spread 
number. In the finite prime field ¥ p of p elements, there arc (p + 3) /2 or (p+l) /2 spread numbers, depending 
respectively on whether p is congruent to 3 mod 4 or 1 mod 4. 

1.3 Definition of the spread polynomials 

The spread polynomials S n (s) arise when we consider the rational analog of multiples of an angle. They may 
be defined recursively, over a general field, without any reference to geometry as follows: 

50 (s) = 

51 (s) = s 

S n (s) =2(1- 2s) S n _! (s) - S n _ 2 (s) + 2s. (1) 

The coefficient of s n in S n (s) is a power of four, so over any field not of characteristic two, S n (s) is a polynomial 
of degree n. The generating function for the spread polynomials was computed by M. Hirschhorn, it is 

S(s,t) = yS n (s)t n = - w t5(1 + t) , -. (2) 

^ (l-t)(l-2t + t 2 + 4ts) w 

The spread polynomials arc intimately linked to geometry, in the following sense. If two intersecting lines Iq 
and l\ in the plane make a spread of s, then we may reflect Iq in l\ to obtain l 2l reflect l\ in l 2 to obtain I3, and 
so on. We will see that the spread s (lo, l n ) is then S n (s). For this description we only need the lines through the 
origin, so it is really a statement about the associated projective line. In fact our statement of this fact, the Spread 
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of a power theorem, will be stated in sufficient generality to include also some non-Euclidean geometries, and will 
involve a multiplicative structure on the non-null points of the projective line. 

Another key fact about the spread polynomials, which connects with their geometric interpretation, is the 
following. 

Theorem (Spread composition) For any natural numbers n and m, 

We will give a proof later using the relations between spreads and rotations; for another see page 110]. 

1.4 Table and graphs of spread polynomials 

Here are the first few spread polynomials. Note that S2 (s) is the logistic map. 



So 


s) 


= 




Si 


I s ) 


= s 




62 


s) 


= 4s - 4s 2 = 4s (1 - s) 




s 3 


's) 


= 9s - 24s 2 + 16s 3 = s (3 - 4s) 2 




S4 


s) 


= 16s - 80s 2 + 128s 3 - 64s 4 = 16s (1 - 


s) (1-2 


s 5 


s) 


= 25s - 200s 2 + 560s 3 - 640s 4 + 256s 5 


= s (5 - 


S6 


I s ) 


= 4s(l-s) (3-4s) 2 (l-4s) 2 




s 7 


I s ) 


= s (7- 56s + 112s 2 - 64s 3 ) 2 




Ss 


I s ) 


= 64s (1 - s) (1 - 2s) 2 (1 - 8s + 8s 2 ) 2 




S9 


I s ) 


= s (3 - 4s) 2 (3 - 36s + 96s 2 - 64s 3 ) 2 




S10 


I s ) 


= 4s (1 - s) (l - 12s + 16s 2 ) 2 (5 - 20s - 


f 16s 2 ) 2 



16s 2 



A remarkable fact, already suggested by this list, is that the spread polynomials factor in an interesting way, 
indeed in a more pleasant fashion than the Chebyshev polynomials. We will establish this in a future paper. 

The next Figures show the first eight and twenty five spread polynomials over the real numbers in the range 
< s < 1. Note the interesting ghost patterns that begin to appear as we increase the number of polynomials 
shown; these are related to Lissajous curves, and such a phenomenon occurs also for Chebyshev polynomials. 

Observe also that the spread polynomials arc positive in this range, so do not form an orthogonal family of 
polynomials in the usual sense unless they are translated vertically 



Sn(s) 




Figure 2: The first eight spread polynomials Sq (s) to S7 (s) 
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Figure 3: The first twenty five spread polynomials 

Example 1 Here are some values of spread polynomials in the field F5. The spread numbers in this field are 0, 1 
and 3. 
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S 6 (s) 











1 





Si (a) 
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S 7 (a) 
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3 


4 


S 2 (s) 








2 


1 


2 


Ss(a) 








2 





2 


Ss(a) 





1 





3 


f 


S 9 (a) 





f 





3 


1 


Si(s) 








2 





2 


S10 (s) 








2 


f 


2 


S 5 (a) 





1 


2 


3 


4 


S n {a) 
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2 


3 


4 



The pattern repeats with period 12, that is S n (s) = S' n +i2 (s) for all n and for all s. However for values of s which are 
spread numbers, S n (s) = S n +i (s) for all n, while for values of s which are non-spread numbers S n (s) = S n +e (s) 
for all n. o 

Example 2 Here are the values of the spread polynomials in the field F7. The spread numbers in this field are 
0, 1,3,4 and 5. 



s 
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2 


3 


4 


5 


6 


s 
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3 
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5 


6 


So (a) 























Sxa (a) 
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4 


3 


6 


Si (s) 
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4 


5 


G 


Su (s) 
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4 


1 


4 


G 


S 2 (s) 








G 


4 


1 


4 


G 


Sis {a) 





1 


1 


3 


4 


5 





S3 (a) 





1 


1 


■ r j 


4 


3 





Sw (s) 








G 











(i 


Si(a) 








G 


1 





1 


G 


5l7 (*) 
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5 


G 


S 5 (a) 





1 


2 


5 


4 


3 


G 


Sis (a) 











4 


1 


4 





Se(a) 











4 


1 


4 





S19 (a) 





1 


2 


5 


4 


3 


G 


S 7 (a) 





1 


2 


3 


4 


5 


G 


S20 (a) 








G 


1 





1 


G 


S s (a) 








G 











G 


S21 (a) 





1 


1 


5 


4 


3 





S 9 (a) 





1 


1 


3 


4 


5 





S22 (a) 








G 


4 


1 


4 


G 


S10 (s) 








G 


4 


1 


4 


G 


S 23 (s) 





1 


2 


3 


4 


5 


(i 


S n {a) 





1 


2 


5 


4 


3 


G 


#24 (a) 























S12 (a) 











1 





1 





S25 (s) 
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2 


3 


4 


5 


G 



The pattern repeats with period 24, that is S n (s) = 5 n +24 (s) for all n and for all s. For values of s which are spread 
numbers, however, S n (s) = S n +s (s) for all n, while for the non-spread numbers S n (s) = S n +e (a) ■ o 

1.5 Spread polynomials evaluated at s = 1/3 

Consider the spread s = 1/3 formed by the long side Iq and the diagonal h of an AA sheet of paper, as discussed 
previously. The numbers S„ (1/3) = s (Iq, l n ) are the spreads formed by successive reflections, or rotations, of these 
two lines. 
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The recursive formula |T]) allows us to calculate the following list of prime power factorizations: 



Sr (1/3) = 


= 3- 1 


S 2 (1/3) = 


= 2 3 3 -2 


S 3 (1/3) = 


= 5 2 3~ 3 


5 4 (1/3) = 


= 2 5 3~ 4 


S 5 (1/3) = 


= 3~ 5 


S 6 (1/3) = 


= 2 3 5 2 3~ 6 


S 7 (1/3) = 


= 43 2 3 -7 


S 8 (1/3) = 


= 2 7 7 2 3 -8 


S 9 (1/3) = 


= 5 2 19 2 3 -9 


S w (1/3) 


= 2 3 11 2 3 -10 


Su (1/3) 


= 197 2 3 -11 


S12 (1/3) 


= 2 5 5 2 23 2 3~ 12 


Sis (1/3) 


= 1249 2 3 -13 


Su (1/3) 


= 2 3 13 2 43 2 3~ 14 


Sis (1/3) 


= 5 4 29 2 3 15 


5i 6 (1/3) 


= 2 7 17 3 


Sir (1/3) 


= 9791 2 3- 17 


Sis (1/3) 


= 2 3 5 2 19 2 73 2 3- 18 


Si9 (1/3) 


= 26107 2 3~ 19 


^20 (1/3) 


= 2 5 11 2 241 2 3- 20 


S21 (1/3) 


= 5 2 43 2 167 2 3- 21 


S 22 (1/3) 


= 2 3 197 2 263 2 3- 22 


S 23 (1/3) 


= 139 2 2207 2 3~ 23 


5 24 (1/3) 


- 2 7 5 2 7 2 23 2 47 2 3- 24 


S 25 (1/3) 


= 149 2 1949 2 3~ 25 


S 26 (1/3) 


= 2 3 131 2 1249 2 3~ 26 


S27 (1/3) 


= 5 2 19 2 53 2 433 2 3~ 27 


S 2S (1/3) 


= 2 5 13 2 43 2 1511 2 3- 2 


^29 (1/3) 


= 6973919 2 3" 29 


S 30 (1/3) 


= 2 3 5 4 11 2 29 2 239 2 3" 



Let's make some empirical observations about the above table. For each n, S n (1/3) is a fraction whose denom- 
inator is 3" . The numerator is divisible by 2 precisely when n is even, and the power of 2 appearing is odd. The 
other factor of the numerator is a square. Prime factors of the numerator seem to occur periodically. For an odd 
prime p 7^ 3, define m (p) to be the smallest natural number m such that for all positive multiples n of m, S n (1/3) 
has a factor of p, if such an m exists. From the table, we may guess that this number for small primes p is: 

771 (5) — 3 to (7) = 8 m(ll) = 10 to (13) = 14 
to (17) = 16 to (19) =9 to (23) = 12 to (29) = 15. 

However larger primes also appear in the table, for example perhaps 

to (6973919) = 29. 

The aim of this paper is to try to begin to explain these numbers, and to show that the phenomenon is not dependent 
on the initial spread s = 1/3. 

We adopt the convention that a rational number a is divisible by a prime p precisely when a = pj3 with (3 a 
rational number which can be expressed as (3 = c/d with c and d integers, and d not divisible by p. In this case we 
say p is a factor of a. The following is the main result of this paper. 

Theorem (Spread Periodicity) For any rational number s = a/b and any prime p not dividing b, there is 
a natural number m such that S n (s) is divisible by p precisely when m divides n. This number to is a divisor of 
either p — 1 or p + 1. 

Corollary For any rational number s = a/b, any prime p not dividing b occurs infinitely often as a factor of 
the numbers S n (s) for n = 1,2, 3, ■ ■ ■ . 

In addition, if we find a prime p appearing as a factor of Sk (s), then we can be sure that p will appear as a 
factor of any spread S n k (s), for n = 1, 2, 3, • ■ ■ . For example, from the above observations it follows that (1/3) 
is divisible by 6973919. Note however that we arc not able to address the more difficult problem of determining 
777 (p) for a given p. 

To prove the Spread periodicity theorem, we will explore the metrical geometry of the one- dimensional projective 
line, and show that finite geometries, both Euclidean and non- Euclidean, play a role. 

2 Geometry of the projective line 

The metrical structure of one-dimensional geometry over a general field F, not of characteristic two, was investigated 
recently in [7j. There are two distinctly different contexts: affine and projective, and in this paper it is the projective 
setting that is of primary interest, because we arc interested in rotations and reflections of one-dimensional subspaces 
of a two-dimensional space F 2 , which essentially takes place in the one-dimensional projective line. 
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A vector in F 2 will typically be denoted U or V. If U is a non-zero vector, then u = [U] represents the 
corresponding projective point, or p-point for short, namely the line OU through U and the origin O = [0,0]. 
If U = [x,y] then we write u = [U] = [x : y], with the usual convention for proportions that [x\ : y\] = [x 2 : y 2 ] 
precisely when X\y 2 ~ x 2 yi = 0. The projective points constitute the projective line P 1 . 

Now fix a symmetric bilinear form on F 2 , 

[xi,yi] ■ [x 2 ,y 2 ] = axix 2 + b (x x y 2 + x 2 yi) + cyxy 2 (3) 
which is non-degenerate, that is 

ac-b 2 ^ 0. 

The projective point u = [U] is then null precisely when U • U = 0. It is important to realize that such a bilinear 
form on F 2 can naturally provide a metrical structure on the associated projective line P 1 . 

The projective quadrance, or p-quadrance for short, between the non-null p-points u\ = [Ui] and u 2 = [U 2 ] 
is the number 

, (Ui ■ U2) 2 

This is well-defined, and if u\ = [x\ : y{\ and u 2 = [x 2 : y 2 ], then a generalization of the well-known Fibonacci's 
identity 

(xix 2 + yiy 2 ) 2 + (x x y 2 - x 2 y x ) 2 = (x\ + y\) [x\ + y 2 ) 

gives 

{axix 2 + b {x x y 2 + x 2 yi) + cy x y 2 ) 2 



q (u\, u 2 ) = 1 — . 

(axi + 2bx\y\ + cyf) (ax\ + 2bx 2 y 2 + cy^) 

(ac - b 2 ) (xiy 2 - x 2 yi) 2 



(ax\ + 1bx\y\ + cyf) {ax\ + 2bx 2 y 2 + cy%) ' 

This shows that 5(1*1,1*2) = precisely when iti = 1*2. The pair (P 1 ,^) denotes the projective line together with 
the projective quadrance q. 

The following is the fundamental formula for projective trigonometry in such a one-dimensional setting, with a 
proof similar to the one in [5] . In planar rational trigonometry this law is called the Triple spread formula, and is 
the analog of the fact that the sum of the angles in a triangle in the real Cartesian plane is ir. 

Theorem 3 (Projective triple quad formula) If u = [U], v = [V] and w = [W] are non-null p-points, then 
the p-quadrances q w = q (u, v), q u = q (v, w) and q v = q (it, w) satisfy 

(q u + 1v + q w ) 2 = 2 (ql + q 2 v + q 2 w ) + iq u q v q w . 

Proof. If we write ajj = U ■ U and bjjv = U ■ V then 

_ ayaw — by W _ ajjaw — b\j W _ a\jay — b\j V 

q-a — Qv — ana q w — . 

ayayy au a w ajjay 

The following is an algebraic identity in the abstract variables an, ay, aw, bjjv, buw and by\y- 

{q u + q v + qwf - 2 (ql + ql + ql) - 4q u q v q w 

= (ajjayaw + 2buvbuwbyw — a ubyw ~ o.yblj W — awb uv ) 

x (2buybuwbyw — auayaw + auby W + aybfr W + awbfjy) a^ a v 2 a u 2 . 

But the first factor on the right hand side is the determinant 



ay bjjv bjjw 

bjjy ay by\y 

buw byw aw 



U-U U -V U-W 
U ■ V V ■ V V-W 
U-W V-W W -W 



which is zero since U, V and W arc coplanar. 
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Example 4 For an integer k, the symmetric bilinear form 

(xi, y\) ■ (x 2 , y 2 ) = xix 2 + kyiy 2 

is non- degenerate provided that the characteristic o/F does not divide k, which we henceforth assume. In particular 
we assume that k is non-zero. If u\ = [x± : y±] and u 2 = [x 2 ■ y 2 \ then the associated p-quadrance qk is 

1 v 1 (a;ix 2 + ky x y 2 f 

q k {ui,u 2 ) = 1 - 2 2 2 2 



(arf + fcyf ) + fcy 2 2 ) 
k (xiy 2 - x 2 yi) 2 

~ (a?? + fcy?) (^1 + feyl) ' 
3 Isometries of the projective line 

Let's now show how spread polynomials link to reflections and rotations in one-dimensional geometry. Suppose q 
is some fixed choice of p-quadrance on P 1 . An isometry of (P 1 , q) is a map r : u — > ut that inputs and outputs 
non-null p-points, and satisfies for any non-null p-points u and v 



For a 2 x 2 matrix 



q (u, v) = q (ut, vt) . 



representing the following linear transformation on F 2 : 

a b 



[x, y] ->■ [x, y] 



c d 



[ax + cy,bx + dy] 



define the corresponding projective transformation r 

[x : y] t = [ax + cy : bx + dy] 

and denote it by 



a b 
c d 



The matrix for such a projective transformation is determined only up to a scalar, so that 



a b 




Xa 


\b 


c d 




Ac 


Xd 



for any non-zero number A. We adapt the following from [8]. 

Theorem 5 (Isometries of the projective line) An isometry of the projective space (P 1 , qk) must be a projec- 
tive transformation of one of the following types: 



P[a:b] 



for some non-null projective point [a : b] . 



a b 
-kb a 



\a:b] 



a b 
kb —a 



Proof. Suppose that r is an isometry of (P 1 , qk) and that it sends i\ = [1 : 0] to [a : b] and i 2 = [0 : 1] to [c : d]. 
Since i\ and i 2 are non-null p-points, [a : b] and [c : d] must also be non-null. Since qk (ii,i 2 ) = 0, we must have 

ac + kbd = 0. 

So [c : d] = [kb : —a] . Now given an arbitrary non- null projective point u=[x:y] with ut = v = [w : z], 



qk (u,ii) 



ky 2 



x 2 + ky 2 



= q k (v, [a : b}) 



k [az — bwY 



(a 2 + kb 2 ) {w 2 + kz 2 ) 
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and 



qk (u,i 2 ) 



x 2 + ky 2 



q k (v, [kb : -a}) 



(aw + kbz) 



(a 2 + kb 2 ) (w 2 +kz 2 )' 



Comparing these two equations gives 



so that either 



x : y = (aw + kbz) : (az — bw) 



x : y = (aiv + kbz) : (az — bw) or x : y = — (aw + kbz) : (az — bw) 
Using matrix notation, either 



[x:y] = [w: z] 



a —b 
kb a 



or [x : y] — [w : z] 



-a —b 
-kb a 



a b 
—kb a 



a b 

kb —a 



Inverting gives either 
If [a : b] is non-null, then both of these are isometrics, due to the equations 

and 



a b 




"1 


0" 




a b 


T r 


—kb a 







k 




—kb a 





a 2 + b 2 k 




1 


0" 


a 2 k + b 2 k 2 







k 



a 


b ' 




"1 


0" 




kb 


—a 







k 





We call 



P[a:b] 



a b 
kb —a 

a b 
—kb a 



'a 2 + b 2 k 
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and cr [a:fc ] = 
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b ' 






kb 


—a 







respectively a projective rotation and a projective reflection. These realizations as projective matrices extend 
isometries also to null points. 

The identity transformation 

"l 0" 



P[l:0] = 

is a projective rotation. 

Our convention for compositions is u (t\T2) = (uti) t%. 



1 



Theorem 6 (Composition of isometries) For any non-null p-points [a : b] and [c : d] 



°~[a:b}°~[c:d} ~ P[ac+kbd:ad-bc] P[a:b]P[c:d] ~ P[ac-kbd:ad+bc] 

P[a:b}^[c:d\ = 0~[ac+kbd:ad-bc] °~[a:b}P[c:d] = < J [ac-kbd:ad+bc}. 

Proof. This is a straightforward verification. Another generalization of the Fibonacci identities, 
(ac + kbdj 2 + k(ad- be) 2 = (a 2 + kb 2 ) (c 2 + kd 2 ) = (ac - kbd) 2 + k (ad + be) 2 

shows that the resultant isometries are also associated to non-null points. ■ 

Define G = G (qk) to be the group of isometries of qk, and distinguish the subgroup G e = G e (qk) of projective 
rotations pu-u ■ These latter are naturally in bijection with the non-null p-points. From the Composition of isometries 
theorem, the subgroup G e is commutative. 

The cosct G = G (qk) consists of projective reflections <7r a: M, and these too are also naturally in bijection with 
the non-null p-points. The group G naturally acts on the space of non-null p-points. This action is transitive since 
crr a . 6 ] and P[ a -b] both send [1 : 0] to [a : b]. 

Since projective rotations are in bijection with non-null p-points, we can transfer the multiplicative structure 



P[a:b]P[c:d] — P[ac~kbd:ad+bc] 
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of projective rotations to non-null p-points. So for it = [a : b] and v = [c : d] non-null p-points, define their fc- 
product by the rule: 

uv = [a : b] [c : d] = [ac — fc6d : ad + be] . 

When k = 1 this multiplication is familiar from the two-dimensional setting of complex numbers. Note however 
that we are here working in the one-dimensional situation, over a general field, and allowing different values of k. 
The resulting group is commutative, has identity 

e = [1 : 0] 

and the inverse of u = [a : b] is it -1 = [a : —b}. For a non-null projective point u, we let it™ = uu ■ ■ ■ it (n times) 
denote the n-th power of u. Of course this depends on the prior choice q~qkoi p-quadrance. 



4 Spreads of p-points 

Working in (P 1 ,^), define the /c-spread s k (it) of the non-null p-point u = [a : b] to be the number 

kb 2 



Sk (it) = qk (e, u) 



a 2 + kb 2 ' 



Then Sk (it) = precisely when u = [1 : 0] = e, and (u) = 1 precisely when it = [0 : 1]. Since multiplication by v 
is the same as applying the projective rotation p v , multiplication is an isometry, so that for any non-null p-points 
v and it, 

kb 2 

The next result connects spread polynomials and powers. 

Theorem 7 (Spread of a power) In (P 1 ,^), if u is a non-null projective point with Sk (it) = s, then for any 
natural number n, 

Sk (it™) = S n (s) . 

Proof. Suppose that it = [a : b] so that 

kb 2 

Sk (u) = ■ 



2 + b 2 k 

We know from the Isometries of the projective line theorem that the projective rotation p u has the form 



Pu 



a b 
-kb a 



We diagonalize the matrix ( * ^ J using a number r satisfying r l = —k. If the field F does not contain such a 



,2 

number, the following equations take place in the quadratic extension F (r). First verify that 

-l 



a b\ _ f 1 — 1\ fa + br Ul -1 
a) \r r I \ a — br I \r r 



Then for any natural number n, 



-l 



a b\ = 1 -1\ (a + br) n A -1 N 

-kb a J \r r J \ (a — br) n J \ r r 

\ (a + brf + I (a - br) n i (a + 6r) n - 4 (a - 6r)' 
v -| (a + brf - | (a - 6r) n | (a + br) n + | (a - 6r) n 

Taking the proportion determined by the first row gives 

u n = [r (a + br) n + r (a — br) n : (a + br) n — (a — br) n ] 
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t n = s k (u n ) 

where A = a + br and B = a — br. 
Clearly t = 0, and since r 2 = — k, 



((a + br) n - (a - br) n f _ (A n — B n ) 2 
4 (a + br) n (a - br) n ~ 4A n B n 



kb 2 

U 



a 2 + b 2 k 

To show that t n = S n (s) for all natural numbers n, we establish the identity 

£ n -2(l-2s)t n _i+t n _ 2 -2s = (4) 
for all n > 2. The left hand side of flU is 

(4" — B n ) 2 ( 2{A-B) 2 \ (A"- 1 -^"- 1 ) 2 (A"' 2 -B n - 2 ) 2 2 (A - B) 2 
AA n B n + y + 4AB J AA^B™- 1 4^-2 B n-2 + J^B 

and if we remove a factor of (4A n B n )~ 1 this becomes 

- A 2n + 2A n B n - B 2n + (A 2 + B 2 ) ^ 2( ™- 1} - 2A n - 1 B' 1 - 1 + B 2 ("- 1 )^) 

- A 2 B 2 (A 2{n - 2) - 2A n - 2 B n - 2 + B 2{n - 2) ^j + 2A n - 1 B n - 1 (A 2 - 2AB + B 2 ) , 
which after expansion is identically zero, so holds independent of the particular choices of A and B. ■ 
Theorem 8 (Spread composition) For any natural numbers n and m, 

S n o S m , — S n7n . 

Proof. Working over the rational numbers with the Euclidean quadrance q, the Spread of a power theorem shows 
that if s = s (u) for some p-point u, then s (u m ) = S m (s) . It follows that 

S n (S m (*)) = s ((u m ) n ) = s (u mn ) = S nm (s) . 

Since this holds for more than nm different values of s, and both S n o S m and S nm are polynomials of degree nm, 
we conclude that 

5 Spread numbers in a prime field 

Fix an integer k which is not zero in the field ¥ p of characteristic p ^= 2, in other words which is not divisible by 
p. Define a number a in ¥ p to be a fc-spread number precisely when a (1 — a) is k times a square. If k = 1 this 
agrees with our earlier usage. 

For p an odd prime, the finite prime field F p contains an equal number of non-zero squares and non-squares. 
This follows for example from the standard fact that the multiplicative group of a finite field is cyclic. Note that 
for p = 3 mod 4 every number is either a square or the negative of a square, since —1 is not a square, and so every 
number is either a 1-spread number or a (— l)-spread number. However for p = 1 mod 4 the 1-square numbers and 
the (— l)-square numbers agree. 

Every number a in a field is a fc-spread number for at least one k, for example and 1 are /c-spread numbers 
for all fc, and for other a we may choose k = a (1 — a). The next result generalizes the Spread number theorem in 
d Theorem 34]. 

Theorem 9 (Spread number) For any non-null p-points u and v in (P 1 ,^), the p-quadrance q k (u,v) is a k- 
spread number, and conversely for every k-spread number q there exist non-null p-points u and v with q k (u, v) = q, 
and so there exists a non-null p-point w with s k (w) = q. 
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Proof. If u = [xi : yi] and v = [x 2 ■ J/2] then 



_> 



(x( + ky() (xj. + kyi) 
in which case 

M \ , (xix 2 + ky^f (x 2 y! - x^f 

q (1 - q) = K 2 2 

[x\ + ky 2 ) {x\ + kyf) 

which is k times a square, so q is a /c-spread number. Conversely suppose that q is a fc-spread number, so that 
g(l — q) = kr 2 for some number r in the field. If q = 1 then it is the p-quadrance between [1 : 0] and [0 : 1]. 
Otherwise the p-quadrancc between the p-points u = [1 : 0] and v = [1 — q : r] is 

r 2 q (1 - q) 

— q. 



(1-qY + kr 2 (l-qy + q(l-q) 

Note that 

(1 - qf + kr 2 = (1 - qf + q (1 - q) = 1 - q 

is indeed non-zero so both u and v are non-null. 

Once we have u and v we can multiply both by u~ 1 to obtain 1 and w 

q k (l,w) = s k (w) = q. ■ 

To determine fc-spread numbers in a prime field ¥ p it suffices to know the squares in the field. In terms of the 
Legendre symbol 

f a\ J 1 if a is a square in ¥ p 

\p J \ — 1 otherwise 

<j is a fc-spread number precisely when 



p J \P J xP) V p 

So it suffices to know whether or not q, q — 1 and — k are squares. If k = 1 (Euclidean geometry) then the above 
gives a straightforward recipe for finding spread numbers from square numbers. 

Example 10 In F13, ( ^\ = 1 cmc? since the square numbers are 0, 1,3,4,9, 10 and 12, the spread numbers are 
q = 0, 1, 4, 6, 7, 8 and 10, as t/iese are t/iose g f aside /rom and 1 ) whose Legendre symbol agrees with that of q — 1. 
Whereas in Fn, ^ = —1 and since the squares are 0, 1, 3, 4, 5 and 9, the spread numbers are 0, 1, 2, 3, 6, 9 and 
10, t/iose g (aside from and 1^) whose Legendre symbol disagree with that of q — 1. o 



6 Spread periodicity 

Suppose F = F p for some odd prime p. In (P 1 , g^) there are exactly p+ 1 p-points, but the number of null p-points 
depends also on fc; there are null p-points precisely when 

x 2 +ky 2 = 

has non-zero solutions. This occurs precisely when — k is a square modulo p 1 and in this case there are exactly 2 
null p-points, and so the order of the group of rotations G e is p — 1. Otherwise there are no null p-points and G e 
has order p + 1 . 

We regard an element of F p as an expression of the form a/b, where a and b are integers with 6 not divisible by 
p, with the convention that a/b = c/d precisely when 

ad — be is divisible by p. 

Thus every rational number a/b with 6 not divisible by p represents also an element of F p . So the values of the 
spread polynomial S n (s) over ¥ p are just obtained by reducing mod p. 
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Theorem 11 (Spread periodicity) Fix an odd prime p. For any element s in ¥ p , there exists a positive integer 
m with the property that S n (s) = in ¥ p precisely when n is a multiple of to. Furthermore m divides either p — 1 
or p + 1. 

Proof. Choose a non-zero integer k such that s is a fc-spread number in F p . The Spread number theorem then asserts 
that there is a non-null projective point u in the projective line (P 1 , qk) over F p with Sk (u) = s. Multiplication by 
u is a projective rotation p u , and by the Spread of a power theorem s& (u n ) = S n (s). So S n (s) = precisely when 
u n = e = [1, 0]. Since p u belongs to the finite commutative group G e , S n (s) = precisely when n is a multiple of 
the order to of p u in G e . This to divides the order of the group G e , which is cither p — 1 or p + 1. ■ 

Corollary 12 For any rational number s = a/b, any prime p not dividing b occurs infinitely often as a factor of 
the numbers S n (s) for n = 1, 2, 3, ■ • ■ . 

Proof. This is an immediate consequence of the theorem. ■ 

To illustrate the theorem for s = 1/3, we will discuss the numbers m (5), m (7) and m (19). 

Example 13 In F5 the squares are 0, 1 and 4, while the 1-spread numbers are 0, 1 and 3. In this field s = 1/3 = 2 
and so s (1 — s) = 3, which is not a square, but s (1 — s) = k x r 2 for k = 2 and r = 2. That means that s is 
a 2-spread number (and hence also a 3-spread number). There are no null p-points for either of the projective 
quadrances qi or q%, since —2 and —3 are not squares in F5. Hence G e has 6 elements in either case, and so the 
order m (5) divides 6. In fact m (5) = 3. A projective point with 2-spread equal to 2 is u = [1 — s : r] = [2 : 1] . Then 
u 2 = [1 : 2] and u 3 = [1 : 0]. o 

Example 14 In Fj the squares are 0,1,2 and 4, while the 1-spread numbers are 0,1,3,4 and 5. In this field 
s = 1/3 = 5 is a spread number since s (1 — s) = k x r 2 with k = 1 and r = 1. Then (P 1 , gi) has no null p-points 
since —1 is not a square, so G e has order 8. Thus m (7) divides 8. In fact m (7) = 8. A projective point with 1-spread 
equal to 5 is u = [1 — s : r] = [3 : 1]. Then the powers of u are u 2 ~ [6 : 1], u 3 = [5 : 1], u 4 = [0 : 1], u 5 = [2 : 1], 
u 6 = [1 : 1] , u 7 = [4 : 1] and u 8 = [1 : 0]. o 

Example 15 In F19 the squares are 

0,1,4,5,6,7,9,11,16,17 

while the 1-spread numbers are 

0,1,2,4,8,9,10,11,12,16,18. 

In this field s — 13 is not a 1-spread number but s (1 — s) = 15 = k x r 2 with k = —1 and r = 2, so s is a (-l)-spread 
number. Then (P 1 ,^-!) has null p-points since 1 is a square, so G e has order 18. Thus m(19) divides 18. In fact 
to (19) = 9. A projective point with {~l)-spread equal to 13 is u = [1 — s : r] = [13 : 1]. Then the powers of u are 
u 2 = [8 : 1], u 3 = [5 : 1], u 4 = [10 : 1], u 5 = [9 : 1], u 6 = [14 : 1], u 7 = [11 : 1], u 8 = [6 : 1] and u 9 = [1 : 0]. o 



6.1 The special spreads 

The spreads 

„ 1 1 3 
s = 0, 1 

' 4' 2' 4' 

play a special role in rational trigonometry. Over the real numbers, they correspond to angles of 0, 30°, 45°, 60° and 
90° respectively, and also to the supplements of these angles. Over a fixed field, any spread polynomial evaluated 
at any one of these five values of s will also yield one of these five values. 

For example if s = 1/4 then S 2 (s) = 3/4, S 3 (s) = 1, ^4 (s) = 3/4, S 5 (s) = 1/4 and S 6 (s) = 0. It is not hard 
to check that 

s 6 0) = 

for any one of these values of s. Since this is true over the rationals, it is true over any other field not of characteristic 
two. 

Example 16 In F7, s = 1/4 = 2 which is a 3-spread number since s (1 — s) = 5 = k x r 2 with k — 3 and r = 2. 
Then 93 has null p-points since —3 is a square, so G e has order 6. Thus to (7) divides 6. In fact to (7) = 6. A 
projective point with 3-spread equal to 2 is u = [1 — s : r] = [3, 1]. Then the powers of u are u 2 — [1 : 1], u 3 = [0 : 1], 
u 4 = [6 : 1], u 5 = [4 : 1] and u 6 = [1 : 0]. o 
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Example 17 In Fn, s = 1/4 = 3 which is a 1-spread number and s (1 — s) = 5 = 4 2 . Then (P 1 , gi) /ias no null 
p-points since — 1 is not a square, so G e has order 12. Thus m (11) divides 12. In fact m (11) = 6. A projective point 
with 1-spread equal to 3 is u = [1 — s : ?*] = [—2 : 4] = [5:1]. Then the powers of u are u 2 = [9 : 1], u 3 = [0 : 1], 
u A = [2 : 1], u 5 = [6 : 1] cmd it 6 = [1 : 0], o 

7 Conclusion 

The results of this paper indicate that there are further rich number theoretical aspects of the spread polynomials. 
Some of these will be studied in future work, particularly the parallels between Spread polynomials and Chebyshev 
polynomials, and the close affinity with cyclotomic polynomials. 
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